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DNA and other biopolymers differ from classical polymers due to their torsional stiff- 
ness. This property changes the statistical character of their conformations under tension 
from a classical random walk to a problem we call the "torsional directed walk". Moti- 
vated by a recent experiment on single lambda-DNA molecules [Strick et al. Science 271 
(1996) 1835], we formulate the torsional directed walk problem and solve it analytically 
in the appropriate force regime. Our technique affords a direct physical determination of 
the microscopic twist stiffness C and twist-stretch coupling D relevant for DNA function- 
ality. The theory quantitatively fits existing experimental data for relative extension as a 
function of overtwist over a wide range of applied force; fitting to the experimental data 
yields the numerical values C = 120 nm and D = 50 nm. Future experiments will refine 
these values. We also predict that the phenomenon of reduction of effective twist stiffness 
by bend fluctuations should be testable in future single-molecule experiments, and we give 
its analytic form. 
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1. Introduction and Summary 

The theory of random walks is one of the most fundamental problems in statistical 
mechanics, with applications throughout physics, biology, and even finance. The discovery 
that polymer conformations afford a concrete realization of this mathematical problem, and 
the understanding that rubber elasticity is inherently an entropic phenomenon, marked 
the birth of polymer physics [[[J. Remarkably, it has recently become possible to apply 
minuscule forces to single molecules of DNA in solution and observe their extension 0. 
Besides allowing a detailed confirmation of the directed random walk model of entropic 
elasticity, these experiments allow direct physical measurement of microscopic (nanometer- 
scale) mechanical properties of DNA relevant to its function, using mesoscopic (micron- 
scale) apparatus. Two linear elastic parameters of DNA have now been measured in this 
way: the bend persistence length A and the intrinsic-stretch modulus 7 [[§ @ H || • 

DNA and other stiff biopolymers differ from classical polymers, however, in that 
they exhibit torsional as well as bend stiffness. Thus their conformations reflect not a 
classical directed walk but a new fundamental problem: the "torsional directed walk" 
(TDW), whose random variables are the direction of each step relative to its predecessor, 
together with a relative axial twist. In this letter we will formulate the version of the TDW 
appropriate to DNA, solve it analytically in a regime appropriate to a recent experiment 
, and show that the model quantitatively fits the data over a wide range of applied forces 
(Fig. 1) !§. Besides being transparent, analytic formulae facilitate systematic least-squares 
fitting to experimental data. We fit to obtain three microscopic elastic constants: the 
twist persistence length C, bend persistence length A, and intrinsic twist-stretch coupling 
D [ |TU[ JITJ . Since A is known independently we have a check on the model. The experiment 
is not sensitive to the other allowed linear-elastic constants such as twist-bend coupling 

We find that the existing data @ yield A = 49 nm, C = 120 nm, and D = 50 nm; future 
experiments will refine these values when fit to our formula. Many authors have sought to 
extract the value of C from both cyclization experiments and fluorescence depolarization 
||13|| . A key point of this paper is that the force regime we study is free from some vexing 
physical and mathematical difficulties which have helped make the determination of C 
from these experiments controversial. In particular, we can use a continuum model with 
no need for the short-length cutoff required to make Monte Carlo calculations tractable 
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Fig. 1: Relative extension of lambda-DNA versus applied force / and 
overtwist a. From top to bottom, the curves are at fixed force 8.0, 1.3, 
0.8, 0.6, 0.3, and 0.1 pN. The dots are experimental data from Fig. 3 
of || , excluding values of /, a where the DNA is known to denature by 
strand separation. Open dots are outside the range of validity of the 
phantom chain model and were not used in the fit. The corresponding 
points from Fig. 2 of |§ were also used in the fit (not shown) , for a total 
of 49 points. The lines are our theoretical predictions after fitting to 
A, C, and D (see text). 



We also give a simple analytical prediction for the reduction of effective twist stiffness 
by bend fluctuations. This renormalization may explain why some other determinations 
of C give slightly lower values than ours. While its existence was appreciated long ago 
by Shimada and Yamakawa [|15|, a clear experimental test has hitherto not been possible. 
While this effect is only marginally visible in the extant data, again future experiments 
should be able to test our prediction by checking the dependence of C eff on the applied 
force (see below) . Our value for D is similar to within the large errors to recent estimates 

13 II- 

In a recent preprint Bouchiat and Mezard have independently addressed several over- 



lapping issues fnj. We comment on their approach below. 
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2. Experiment 

In the experiment of Ref. [|8| , DNA was bound to a wall at one end and a magnetic bead 
at the other, with bonds which did not permit free pivoting. Constraining the orientation 
of the bead with an applied magnetic field thus constrained the orientation of the DNA 
strand at its end. Since the bead was too large for the DNA to loop around it, this 
procedure effectively fixed the total Link of the circuit consisting of the DNA plus a fixed 
imaginary closing path. Rotating the applied magnetic field then allowed the authors to 
freeze the Link to any desired value and find the extension Z for various values of the 
applied stretching force /. Alternatively, the force could be held fixed while the Link was 
varied, as in Fig. 1 above. The intensive strain variable describing Link is the relative 
overtwist a = ALk • (3.4nm/L), where L = 16400 nm is the total contour length. 

Subsequent work showed that the DNA denatures by strand separation for high ap- 
plied stresses, roughly / > 0.4 pN and a < —0.01 or a > 0.03 (D. Bensimon, private 
communication). The loss of twist rigidity is also clearly visible in the curves in ||. We 
have omitted such points from Fig. 1. Also, DNA undergoes a dramatic overstretching 
transition at around 60 pN [|J; all the data discussed here concern forces / < 8pN. 



3. Physical Picture 

We will begin with a picture of DNA as a thin cylindrical elastic rod of fixed contour 
length L; below we will discuss corrections reflecting the more detailed architecture of the 
molecule. The conformations of such a rod under an applied tension are controlled by the 
elastic energy functional 

E/k B T = \ J ^((fij 2 + (fi 2 ) 2 ) + C(0 3 - iv ) 2 ] ds - f ■ Z - 2ttt - Lk . (3.1) 

Here s is arc length, O x 2 (s) are bending strains, fi 3 — uj is twist strain, A, C are the 
bend and twist persistence lengths, / = f /k B T, and u = 27r/3.4nm is the unstressed 



1 Our notation is similar to JT!|. Throughout we will neglect sequence dependence. In 



the force regime in question we expect linear rod elasticity to be a good approximation. 
Higher-order bend elasticity effects are expected to be suppressed by powers of the rod 
radius (1 nm), which is much smaller than any other length scale. Indeed a linear-elastic 
model, the "extensible worm-like chain" , describes accurately the extension of torsionally 
unconstrained DNA up to forces greater than those considered here, with an intrinsic 
stretch modulus more than a hundred times greater than the forces of interest to us . 
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molecule's helix density, r is a dimensionless torque variable which we will choose to obtain 
the required overtwist (Lk) = (1 + a)Lu /27r. 

In the absence of thermal fluctuations, a rod under extensional stress remains straight 
as we apply increasing torque to the ends, then becomes linearly unstable at a critical value 



of torque which increases with the applied force \\L7j . The end-to-end distance Z of the 
rod does not change at all for applied overtwist less than the critical value r crit = 2\J Af, 
since the rod remains a straight line of constant contour length. 

Thermal fluctuations change this picture completely. The rod is never straight; every 
Fourier mode of its shape is excited in accordance with the equipartition theorem of sta- 
tistical physics, so the net length Z is always less than L. The applied tension suppresses 



those modes of wavenumber smaller than q = y f/A, with the dominant contribution 
to 1 — (Z/L) coming from q « q Q . Now when we apply torque below the critical value, the 
modes with the same helical sense as the torque get pushed closer to instability and their 
fluctuations are enhanced, shortening the rod. The modes with the opposite helical sense 
are suppressed, and to first order in a there is no effect, consistent with the obvious fact 
that the elastic energy ( |3.1| ) does not break inversion symmetry. We do expect an effect 
to 0(a 2 ); indeed this is the dominant feature of the data in Fig. 1. As the force decreases 
the effect increases, as seen in the increasing curvature of the curves in the figure. 

In eqn. ( |3.1| ) we have neglected any self-avoidance effects; these would appear as 
interactions between rod elements distant in s. In the usual directed walk this is not a 
serious omission: the crossover to self-avoiding-walk scaling occurs only for chain lengths 
much longer, and forces much smaller, than those encountered in DNA. In the TDW 
we must be more careful, since the linking number appearing in is undefined when 
the chain crosses itself. Physically the problem is that the chain can form a loop, pass 
through itself, and in the process lose a unit of Link: the phantom torsional chain cannot 
support any imposed torque. If we seek an equilibrium at nonzero r we must expect to 
find mathematical pathologies; they will arrive in due course. A related problem is that 
the statistical sum for the phantom torsional chain includes all knotted configurations, an 
error with noticeable effects [[18] . 



One approach to this problem is to introduce realistic self-avoidance and knot rejection 
into eqn. ( |3.1| ) |[18|||14|| . The resulting nonlocal model requires numerical Monte Carlo 
solution, and the results depend on the details of the chain interaction chosen. From the 
physical picture, however, it is clear that at high enough applied tension / the problematic 
loops and knots will be so rare as to be negligible: the chain remains nearly straight, and 
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we may use (3J.) without modification. We will find below the precise condition to be in 
this regime; it corresponds to the solid dots in Fig. 1. 

Eqn. ( P?T| ) also neglects any effects of rod anisotropy. For example, bending into 
the major groove ( "Roll" ) is obviously easier than bending in the perpendicular direction 



("Tilt"); less obvious is an allowed twist-bend coupling [^2| . Such anisotropies can lead to 
chiral effects, for example an asymmetry between a and —a, but their effects on Fig. 1 are 
negligible (J.D. Moroz and P. Nelson, in preparation). Indeed we only expect the helical 
pitch to affect entropic elasticity when the dominant wavenumber q approaches w , i.e. 
at unattainably large forcesJl 

Instead the asymmetry visible in Fig. 1 has its origin in the intrinsic stretch elasticity 
of a chiral rod, which gives Z a contribution linear in a and independent of / [^TJ. This 
coupling may be relevant for the binding of the protein RecA to DNA, which stretches and 
untwists the DNA The effect should be masked by the entropic a 2 term for small 

force, emerging when the latter is suppressed at high force, exactly as seen in Fig. 1. 



4. Calculation 

We now sketch a calculation embodying the above physical picture (J.D. Moroz and P. 
Nelson, in preparation). We introduce three local configuration variables: a unit vector t(s) 
describing the tangent to the chain, and an angle ((s) for the remaining torsional degree 
of freedom. We take the applied force along the z direction and orient the chain so that 
t=z in equilibrium. To define ((s) we use Fuller's local formula for the writhe of a curve 
whose tangent never points along the -z axis ||20||||21|| : Wr= J ds(t x ^|) • z/(l + t ■ z). 
Combined with White's theorem that 2?rWr + Jfi 3 ds is a topological invariant, we see 
that fi 3 + (t x 4 5 )-z/(l+t-z) must be a total derivative; we will call this quantity uj + ^ 
and eliminate fi 3 (s) in favor of C( s )-^ The advantage of this choice is that ( |3.1|) is now 
quadratic in £, which may be summarily eliminated. 

2 Certainly the omitted anisotropic couplings will also renormalize the constants A, C 
in ( p.l|) (J.D. Moroz and P. Nelson, in preparation). The model ( |3.1|) is to be regarded as 
coarse-grained to the scale of the helix pitch. 

3 A more elegant, equivalent, approach takes the configuration variables to be a 3 x 3 
rotation matrix; the torque term then takes the form — r J ds[fi 3 + fi 3 ]/(l + t • z), where 
fi^ are the space-fixed angular velocities of a rigid body (J.D. Moroz and P. Nelson, in 
preparation) . 
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It proves convenient to define dimensionless quantities K 
Ks/A; we then find 



Af— r 2 /4 and s 



E/k B T = S + 



ds 



|ii|2 , 2fA 



(1-z-t) 



2t 



(t x dt/ds) -z/(l + t- z) 



(4.1) 



where £ = —L(f + r 2 /2C) and dot denotes The second term defines a nonlinear 
fluctuation problem, which we will expand in powers of 1/K. From its partition function 
Z(f,r) we may then extract the extension and excess link as 



d_ 

Of 



logZ , u (l + a) 



L !h 



log Z 



(4.2) 



We will use the second of these to solve for r(f, a), then substitute into the first to get the 
desired extension Z(f, a). 

To find Z we adapt the standard trick used in the wormlike chain 0[]5|: for long 
chains Z approaches the unnormalized correlation function of t. Holding t(0) fixed, this 
correlator i/j(t,s) obeys the Schrodinger-like equation ip = —Htfj, where (J.D. Moroz and 
P. Nelson, in preparation) 

H = K(l - cos6) - ^L 2 + £(1 + COS0)- 1 L z + - (1 + cos #) -1 ) (1 - cos#) . (4.3) 

Here cos 6 = t ■ z, L 2 is the angular part of the Laplace operator, and L z is the azimuthal 
derivative. The main novelty of this derivation is the presence of first-order derivatives in 
( [Lip when r / 0, leading to the r 2 terms in K and in ( |4.3| ). The asymptotic value of Z 
is then controlled by the lowest ("ground-state") eigenvalue of (|4.3|) via Z oc e _£ - 0_Aoi . 

Unfortunately ( |4.3| ) has no ground state for any nonzero r, due to its singularity at 
6 = 7r! This unphysical pathology was predicted in the previous section; its mathematical 
origin is the breakdown of Fuller's formula when t = — z. To see the connection to the 
physical discussion, note that the unphysical link-dropping process in the phantom tor- 
sional chain necessarily involves the tangent t(s) passing through -z at some intermediate 
point. As discussed above, a physically meaningful and analytically tractable resolution 
to the problem is to restrict attention to large /. We can then solve (|4.3| ) in perturbtion 
theory about 9 = 0, where the problem is invisible, provided the perturbative ground state 
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value is smaller than the "tunneling barrier" of ( |4.3| ) .□ Imposing this condition and K 2 > 2 
selects the solid dots in Fig. 1. 

The perturbative solution of ( |4.3| ) gives the ground state eigenvalue A = 1 — 1/4K — 
1/6AK 2 ---, so O gives a = ^{C' 1 + {AKA)' 1 + 0(K~ 3 )), Z/L = 1 - (2K)~ 1 (1 + 
1/64K 2 + 0(K~ 3 ). Solving by iteration gives the torque 



r(f, a) = u a/ [C' 1 + (4A\J Af/k B T) , (4.4) 

plus corrections of 0(K~ 3 ). Formula ( |4.4| ) displays the promised renormalization of twist 
stiffness by bend fluctuations. While direct torque measurements are not currently possible, 
this effect nevertheless enters the force curves since r enters Z/L. 

Assembling the pieces gives our theoretical prediction for the force curve: the relative 
extension Z/L is 

Z(f, a)/L = 1 - \(Af/k B T - (r/2) 2 - 1/32)- 1 / 2 + (/ - Da)/ 1 + A/LK 2 , (4.5) 

with r from ( |4.4| ). We have improved the previous formula (J.D. Moroz and P. Nelson, 
in preparation) by summing perturbation theory, introducing the intrinsic-stretch terms 
mentioned in the previous section, and including a small finite-length correction. 

In formulas (|4.4|) (|4.5|) the parameters w , L are known and we use the value 7 = 
1100 pN for the linear stretch constant obtained from higher- force experiments 0. This 
leaves A, C, and D which we fit to the experimental data after the cuts described above. 
Fig. 1 shows that a single choice of A : C, D fits all the curvesJi A nontrivial check is that 
several points just outside our accepted set, not used in the fit, nevertheless lie on our 
theoretical curves. 



4 To justify perturbation theory itself we note that it gives an excellent approximation 
to the exact solution of the wormlike chain when K > 1, as may be expected from the 
form of the leading anharmonic correction below. Nonperturbative effects in the variational 
approach to the WLC J5| are also small when K > 1. Note that our data cuts also eliminate 
the region a > 2-k/CuJq where plectoneme formation (and hence large self-avoidance effects) 
is expected . Raising our threshold on K selected fewer points with little effect on our 



result. 

5 The variance of the data from our curves is = 0.013, comparable to the visible 

scatter in the data. The formal covariances for A, C, Q correspond to very small errors; in 
practice the fit is visibly worse for C outside the range 70 < C < 150. 
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Bouchiat and Mezard have recently analyzed the / = 0.1 pN curve [16]. They inde- 
pendently obtained formulae equivalent to (|4.2| ) (fOj). Since the 0.1 pN curve is entirely 
outside the range of validity of the phantom chain model, however, they were obliged to 
address the unphysical pathology of (|4.3|) , both by Monte Carlo simulation and by intro- 
ducing a new short-scale cutoff. The new cutoff introduces a new unknown parameter into 
the theory, and moreover does not correspond in a simple way to the actual physics of 
self-avoidance. Nevertheless they found an impressive fit to the 0.1 pN curve and relative 
insensitivity to their choice of cutoff for |a| < 0.015. 

Formulas (|4.4|) ( f4.5|) make many predictions. We note that the fit value of A is con- 
sistent with jTj. Omitting either the f/^y or the —Da/'y terms makes a poorer fit, as does 
replacing (|4.4|) by the naive r = u a/C eS with constant C eff . (It is interesting to note that 
the value C eS = 98 nm obtained in this way is closer to earlier determinations than our 
microscopic value of C.) Definitive test of Q4.4|) must await further experiments. 
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